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ABSTRACT 

We present a detailed investigation of SBS 1150+599A, a close binary star hosted by the 
planetary nebula PN G135.9+55.9 (TS01, Stasihska et al. 2009). The nebula, located in the 
Galactic halo, is the most oxygen-poor one known to date and is the only one known to harbor a 
double degenerate core. We present XMM-iVewfon observations of this object, which allowed the 
detection of the previously invisible component of the binary core, whose existence was inferred 
so far only from radial velocity and photometric variations. The parameters of the binary system 
were deduced from a wealth of information via three independent routes using the spectral energy 
distribution (from the infrared to X-rays), the light and radial velocity curves, and a detailed 
model atmosphere fitting of the stellar absorption features of the optical/UV component. We find 
that the cool component must have a mass of 0.54 ± 0.2 M Q , an average effective temperature, 
T cff , of 58 000 ± 3 000K, a mean radius of 0.43 ± 0.3 R , a gravity \ogg = 5.0 ± 0.3, and that 
it nearly fills its Roche lobe. Its surface elemental abundances are found to be: 12 + log He/H 
= 10.95 ±0.04 dex, 12 + log C/H = 7.20±0.3dcx, 12 + log N/H < 6.92 and 12 + log O/H 
< 6.80, in overall agreement with the chemical composition of the planetary nebula. The hot 
component has T ff = 160-180 kK, a luminosity of about ~ 10 4 L Q and a radius slightly larger 
than that of a white dwarf. It is probably bloated and heated as a result of intense accretion and 
nuclear burning on its surface in the past. The total mass of the binary system is very close to 
Chandrasekhar limit. This makes TS 01 one of the best type la supernova progenitor candidates. 
We propose two possible scenarios for the evolution of the system up to its present stage. 

Subject headings: stars: binaries: close: individual(SBS1150+599A): atmospheres: AGB and post-AGB: 
symbiotic: evolution: supernovae; (ISM:) planetary nebulae: individual (PN G135.9+55.9, TS01) 

1. Introduction 
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SBS1150+599A was identified as a planetary 
nebula (PN) in Tovmassian et al. (2001) and sub- 
sequently designated as PN G135.9+55.9. More 
recently, we refer to this object as TS 01 (Stasihska 
et al. 2009). The object has unusually few spectral 
lines for a PN and is renown for its extremely low 
oxygen content (Tovmassian et al. 2001; Jacoby 
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ct al. 2002; Pequignot & Tsamis 2005; Stasiriska 
ct al. 2005, 2009). It is located above the Galactic 
plane at a distance of at least a dozen kpc, which 
places it among a handful of known halo PNe. Di- 
rect images obtained on the ground (Richer et al. 
2002; Jacoby et al. 2002), and most recently by 
HST (Napiwotzki et al. 2005; Stasiriska et al. 2009) 
confirm its PN identification. The observed ex- 
pansion velocity of the nebula (Richer et al. 2003) 
is typical of PNe. But another outstanding feature 
of this PN is that it harbors a close binary sys- 
tem (Tovmassian et al. 2004), revealed serendip- 
itously by the displacement of stellar lines with 
respect to nebular lines within a single observing 
night. Since only one component of the binary 
could be observed in the optical and UV, it was 
suggested that the visible component has a tem- 
perature of 110 000-120 000 K. The lower limit is 
the minimum effective temperature needed to pro- 
duce the observed [NeV] nebular emission line, 
while the upper limit was deduced from the slope 
of the continuum (Tovmassian et al. 2004) . There 
was an ambiguity in the determination of the or- 
bital period, although it was clear that the nucleus 
is a close binary with a period less than 4 hours. 
The high temperature, coupled with high log g, 
determined from the profiles of absorption lines, 
led all studies prior to Stasihska et al. (2009) to 
assume that the observed optical/UV component 
was the central star of the planetary nebula, i.e. 
the post-AGB star that lost its envelope and was 
the source of its ionization. Pequignot & Tsamis 
(2005) suggested that, if the ionizing star were 
even hotter, the deduced oxygen abundance could 
be increased to a more common level for oxygen- 
poor PNe. However, a higher temperature would 
have required a higher reddening to match the 
observed continuum slope, and Tovmassian et al. 
(2004) had already used a higher extinction than 
would normally be estimated for the direction of 
SBS 1150+599A in order to justify a temperature 
of 120 000 K . 

Next, we obtained photometric light curves of 
the binary core of SBS1150+599A (Napiwotzki 
ct al. 2005). The orbital period of the sys- 
tem turned out to be 3.92 hr and, to explain 
the double-humped shape of the light curve, we 
were led to invoke a Roche lobe-filling optical/UV 
component. It was observed that the depths of 
the minima in the light curve are uneven, an 



effect known to occur when the visible compo- 
nent is irradiated by a hotter (more energetic) 
source. The orbital dynamics required that this 
invisible component be another compact object 
of at least 0.85 M (Napiwotzki et al. 2005). Ja- 
coby et al. (2002) pointed out the possibility that 
SBS 1150+599A may be associated with the X- 
ray source 1RXS J115327.2+593959. To detect 
the invisible source of irradiation and reveal the 
other component of the close binary, we observed 
it with the XMM-JVewton X-ray observatory. We 
also conducted new optical spectroscopic observa- 
tions of the object with the Gemini-North tele- 
scope to improve our knowledge of radial veloc- 
ities of the optical/UV component of the binary 
and to better fit photospheric line profiles with 
atmospheric models. We also used the publicly 
available HST STIS observations of the object in 
the UV to bridge the optical and X-ray observa- 
tions discussed here. 

The ionization state and chemical composition 
of the planetary nebula are analyzed in a com- 
panion paper (Stasihska et al. 2009), while here 
we present a multifaceted analysis and modeling 
of the binary system. We analyze the history and 
the future of the stellar system in the light of evo- 
lutionary models for close binary stars. 

In Section 2, we present our new observations; 
in Section 3 we determine the physical parameters 
of the binary; in Section 4 we discuss the evolution 
of the object from the early stages, when it was a 
wide system comprised of main sequence stars, to 
the latest stage of a merging of two white dwarfs 
(WD) with possible type la supernova outcome; 
and in Section 5 we summarize our main results. 

2. Observations. 

2.1. Optical observations. 

A spectrum of TS01, with an ample coverage 
of 3800-9200 A is available in the Sloan Digital 
Sky Survey (SDSS 1 ). It was taken on 2002 May 17 
(spSpec-52411-0953-160). We used the newly cal- 
ibrated spectrum that appears in the SDSS Data 
Release 7. This spectrum provides probably the 
best flux calibration, in a perfect agreement with 
HST spectra (see below). 

The stellar absorption lines of the Balmer series 
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and of He II are difficult to detect due to their shal- 
lowness, the faintness of the object (V~ I8 m ), and 
the presence of very intense emission lines from 
the nebula. Only 7 spectra with measurable ab- 
sorption features were available from our CFHT 
observations (Tovmassian et al. 2004). The short 
orbital period and consequent line smearing by the 
long exposures, and the necessity of relatively high 
resolution to effectively disentangle emission lines 
from absorption indicated the need for observa- 
tions with a larger telescope. 

We proposed to observe TS 01 for a total of 16 
hours, covering four orbital periods, with Gemini- 
North telescope. The observations were scheduled 
for service mode in semester 2006A, but only 20% 
were completed. The available observations were 
performed in two sets: on 2006 May 16 UT, 8 
spectra were obtained, and, on 2006 Jun 09, four 
more were added. The weather conditions dur- 
ing the observations were not ideal. The exposure 
times were 700 s, so, in total, only about 3/4 of 
the orbital period was covered with a resolution in 
phase of 5%. The GMOS spectrograph was used 
with the B1200+G5301 grating, leading to effec- 
tive spectral resolution of 1.65 A (FWHM) and a 
spectral coverage of 3800-5000 A. We observed the 
Balmer series from H^ to the highest members (Hg 
) This spectral range also includes the He II 4686 A 
line, detectable in both emission and absorption. 
Auxiliary images (biases, flatfields, arcs) were used 
to reduce the data using the procedures in gemini 
package within IRAF 2 and prescriptions provided 
by Gemini staff and fellow observers 3 . The stan- 
dard star PG 1545+035, observed with the same 
instrumental configuration on 2006 Aug 30 in ap- 
parently better conditions, was used in an attempt 
at flux calibration. However, the result of this cali- 
bration was not satisfactory. Instead, we used the 
spectrum of TS 01 available in the Sloan Digital 
Sky Survey (SDSS) data base to correct the con- 
tinuum. The Gemini spectra were corrected for 
orbital radial velocity shifts using the orbital so- 
lution described below and then co-added. Com- 
bining the 12 radial velocity (RV) shifted spectra 
allowed us to improve the profiles of stellar ab- 
sorption lines and to get rid of nebular emission 

2 Copyright(c) 1986 Association of Universities for Research 

in Astronomy Inc. 
3 http: / /www. astro. caltech.edu/kelle/gmos/gemini_reduction. html 




i i i i I i i i I i i i I i i i I i i i I i 

3800 4000 4200 4400 4600 4800 

Wavelength (A) 



Fig. 1. — The stellar spectrum of TS01 obtained 
at Gemini North (solid red line). The spectrum is 
a combination of 12 individual spectra, observed 
at different orbital phases, and smoothed with a 
13 pixel boxcar. Before combining, the individual 
spectra were corrected for orbital motion. The 
SDSS spectrum, containing nebular lines, is plot- 
ted in yellow. The black dotted line is the sin- 
gle HST spectrum obtained for the same spectral 
range (see the text). 

lines. The resulting spectrum, after 13 pixel box- 
car smoothing, is presented in Fig. 1. 

Previous observations of TS01 are used here 
to analyze the nature of the stellar core. In ad- 
dition to the above mentioned 7 CFHT spectra 
of lower spectral resolution, they include multi- 
band photometric observations, briefly presented 
in Napiwotzki et al. (2005). They were obtained 
on two consecutive nights with the 2.2 m telescope 
at Calar Alto and the BUSCA CCD camera sys- 
tem that allows simultaneous direct imaging in 
four colors. The differential photometry was per- 
formed using comparison stars in the field of view. 
These photometric data were complemented by 
CCD photometry in the V filter obtained with the 
2.1m telescope at the Observatorio Astronomico 
Nacional in the Sierra de San Pedro Martir (OAN 
SPM) on 2004 Apr 09. Additional photometric 
data were provided by the optical monitoring in- 
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strument (OM) on board XMM during the X-ray 
observations. Optical and UV data from direct 
images in the optical range as well as integrated 
flux from spectroscopic observations were also in- 
corporated into the time series. For the time series 
analysis, the photometric data from the different 
wavelengths and bandpasses were normalized to a 
common mean value and combined. 
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Fig. 2. — The source spectrum in all three detec- 
tors, PN, MOS1, and MOS2 (black, red and green 
crosses respectively) are fitted with the same black 
body model (lines of corresponding colors). The 
photons are concentrated in the lowest 0.1- 0.3 
keV range. There is general consistency between 
all three EPIC detectors at energies below 300 eV. 
In the 0.4 - 0.5 keV bin, the counts start to differ, 
but the statistical significance of this is very low 
taking into account low count rates and calibra- 
tion uncertainties. 

2.2. Ultraviolet observations. 

The first UV data for TS 01 were obtained in 
the far-UV using the Far Ultraviolet Spectroscopic 
Explorer satellite (FUSE) . Details of these obser- 
vations and their results are provided in Tovmas- 
sian et al. (2004). 

Later, observations in the near-UV were ob- 
tained with the Hubble Space Telescope (HST) 
(Obs. ID 9466). In 2003 May, HST obtained 
spectra of TS 01 with the Space Telescope Imag- 
ing Spectrograph (STIS) in FUV, NUV, and CCD 
modes to cover the entire UV and optical. Contin- 
uous (uninterrupted by the Earth occultations) ex- 
posures with the G140L and G230L gratings were 



acquired. Five spectra with each grating and with 
individual exposure times of 4675 s and 2850 s, re- 
spectively, were acquired and combined to produce 
the final spectrum. One 600 s exposure was taken 
of the UV-optical spectrum (G430L grating) to 
connect the UV data with the optical. This last 
spectrum, as a result of its short exposure time, 
fails to reveal relatively weak, though important, 
emission lines in the optical UV, but provides a 
decent stellar spectrum that overlaps nicely with 
the NUV and SDSS spectra. This spectrum is also 
plotted in Fig. 1 and is in very good agreement re- 
garding absorption features. The object was also 
observed in 2003 June with High Resolution Cam- 
era (HRC) of the Advanced Camera for Surveys 
(ACS) in the F334N & F658N filters to obtain im- 
ages in the strong nebular lines of [NeV] and H Q , 
respectively. 

The pipeline-reduced STIS spectra were uti- 
lized to extract the stellar continuum. The inte- 
grated fluxes of individual exposures in UV range 
were also summed to produce a light curve in the 
UV range. 

2.3. X-ray observations. 

TS01 was observed with XMM-JVewton (obs- 
ID 0404220101) on 2006 Nov 01-02 (revolution 
1263) in a continuous 27 ks exposure. The X-ray- 
counting EPIC instruments were operated with 
the thin filter in the PN small window mode and 
full window for the MOS detectors. The object 
was too faint for EPIC-RGS detectors. The op- 
tical monitor (OM) instrument on board XMM- 
Newton took 16 images in a B filter, each of 22 
minutes duration. No pile-up was detected in ei- 
ther of the EPIC detectors. Background photon 
flares were detected during only 35-40% of expo- 
sure, mostly towards its end. The 7.5 hour ex- 
posure is just shy of two orbital periods (2 x 3.9 
hours) of the binary system. The data from the 
first orbit were completely free from background 
flaring effects. The observed mean source count 
rates were 0.033±0.002 in the PN, 0.0025±0.0004 
and 0.0053 ± 0.0005 cnts s" 1 , in the MOS1 and 
MOS2, respectively. 

The data were reduced using XMM-SAS (ver- 
sion 9.0). For the MOS detector, the source and 
background photons were extracted from a circu- 
lar aperture and surrounding annulus correspond- 
ingly. For the PN detector, we tried subtraction 
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of background from 2 different circles near the 
source, since the small window did not allow to use 
a large annulus. We found no substantial differ- 
ences in background removal from different areas. 
Events with detection patterns of up to quadruples 
were selected. 

Background-subtracted spectra in the three 
EPIC detectors and a single blackbody model cor- 
responding to each detector are shown in Fig. 2. 
Estimates of the total galactic H I column density 
varies from iV H = 1.8 X 10 20 (Kalberla et al. 2005) 
to N H = 1.53 x 10 20 cm- 2 (Dickey & Lockman 
1990). The X-ray spectral analysis, with the col- 
umn density fixed to the mean value determined 
for the direction of TS 01 {N H = 1.6 x 10 20 c m - 2 ), 
gives a best fit for kT = 17eV (T« 195 000 K). 
Based upon the XSPEC (Arnaud 1996) modeling, 
the 90% confidence region spans kT = 12.7 — 18.0 
eV. The source is extremely soft and emits only 
in the narrow range spanning 100 - 300 eV. This 
range is notorious for its unreliable calibration 
(e.g. Mateos et al. (2009) for the latest evalua- 
tion) and routinely excluded by observers. How- 
ever there is a good agreement between flux in the 
PN and MOS detectors for 0.1-0.3 keV (Fig. 2). 
Although the PN detector suggests the possibility 
that there is emission in excess of the blackbody in 
the 0.4-0.5 keV bin, the excess is not confirmed by 
either of the MOS detectors, making it unlikely 
that it is a real spectral feature. The analysis 
of ROSAT RASS archival data reveals that the 
source was poorly covered and the background is 
uneven, making spectral fitting useless, though it 
does confirm the extreme softness of the source. 

In Fig. 3, the light curve of the source in the 
PN detector is presented extracted in two energy 
bands, 0.1-0.3 keV and 0.3-10 keV. This light curve 
demonstrates that practically all photons from the 
source are emitted in the narrow soft band, that 
the background flaring occurs mostly in the last 
quarter of the 27 ks exposure, and that the flares 
do not affect the soft band, where the source emits, 
but are rather strong in higher energies, confirm- 
ing their nature as background. A similar pic- 
ture emerges from the MOS detectors. However, 
we have chosen a conservative approach and ex- 
cluded all episodes when the count-rate exceeded 
0.1 counts/s from the analysis of the source in all 
three detectors. 

The X-ray light curve of the source in the PN 
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Fig. 3. — The X-ray light curve of the source in the 
PN detector (top panel). The filled symbols rep- 
resent the light curve in 0.1-0.3 keV range, while 
the open symbols (yellow) are the counts in the 
0.3-10.0 keV band. The count rate above 0.3 keV 
is practically zero, except during the final part of 
the exposure, where they are dominated by back- 
ground flaring. In the extreme soft band, 0.1- 
0.3 keV, the count rate is approximately constant 
(0.033 cnts/sec) for the entire duration of the ob- 
servation. The bottom panel displays the optical 
light curve for comparison. Black points are OM 
measurements on board of XMM-Newton at the 
time of the observations whereas the open boxes 
are V light curve obtained at SPM (2.5 years ear- 
lier; shifted in time according to phases). 

detector shows some flickering but no definite pe- 
riodic variability. There is no correlation between 
the power spectra obtained for the three different 
detectors and the optical light curve. 

The OM measurements in the B filter from the 
pipeline reduction were simply transformed into 
a magnitude scale by taking the logarithm and 
shifted to the same average value as the ground- 
based optical differential photometry. In this way, 
they were used to identify the precise orbital phas- 
ing for the X-ray data. 
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Fig. 4. — The power spectrum of all photomet- 
ric data. The strongest, single peaked maximum 
corresponds to the orbital period, at 6.1159 day -1 . 
The first harmonic, at w 12 days" 1 , is also promi- 
nent in the power spectrum because of the double- 
humped nature of the light curve. 

3. Physical parameters of the close binary. 
3.1. Orbital Period. 

The orbital parameters were roughly deter- 
mined with the discovery of the binarity (Tovmas- 
sian et al. 2004). However the precise orbital pe- 
riod can be determined only with the accumula- 
tion of enough data over a sufficiently long period 
of time. At the moment, we have uneven tempo- 
ral coverage spanning about 1000 days, combin- 
ing multicolor photometry from the ground, op- 
tical monitoring on board of XMM-Newton and 
narrow-band photometry from HST in the opti- 
cal range. We also considered the integrated flux 
in narrow bands from UV spectra obtained with 
HST. All of the photometric data were shifted to 
an arbitrary mean magnitude and analyzed for 
the presence of a period using the discrete Fourier 
transform (DFT) method. The power spectrum 
is presented in Fig. 4 with its strongest peak at 
6.1159 cycles/day, the frequency corresponding to 
the orbital period, and its 2/ alias at 12.23 cy- 
cles/day. The resulting ephemeris is 

HJD = 2452760.6756(5) + 0^163508(3) x E, 

where the zero point T corresponds to the 
deeper minimum in the light curve. The light 
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Fig. 5.— The light curve of TS 01 folded with the 
orbital period. The lower panel presents optical 
light curves. From top to to bottom are: BUSCA 
UV and B, Johnson V from SPM, XMM OM 
(black points), and finally BUSCA R and NIR. 
The upper panel presents HST data: H Q (red), 
[Ne v] (magenta), the integrated spectral flux in 
the NUV (dark blue) and FUV (cyan). The hor- 
izontal error bars in the upper panel denote the 
exposure time converted to orbital phase. 

curves of TS01 folded with the estimated 3.924 h 
orbital period are presented in Fig. 5 in different 
bands. In the bottom panel multi color photome- 
try from Calar Alto is plotted combined with V 
band data from OAN SPM and OM-XMM. In 
the top panel the measurements from a variety 
of HST detectors are displayed. Even though the 
HST data in the optical narrow filters F334N & 
F658N include large contributions from the neb- 
ular emission lines, they nevertheless show vari- 
ability of the stellar core with a similar amplitude 
as in the broad-band filters. The UV light curves 
have the same double-humped shape as their op- 
tical counterparts, but the amplitude decreases as 
the wavelength moves further into UV. The far 
UV observations with the G140L grating have ex- 
posures that are almost a quarter of the orbital 
period, so orbital smearing is severe. Degrading 
the optical light curves to a similar time resolu- 
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tion shows that the small amplitude in the far UV 
light curve is the result of smearing rather than an 
actual change in the amplitude of the variation. 

The possible interpretation of the double 
humped light curve was briefly discussed in Napi- 
wotzki et al. (2005). With the X-ray observa- 
tions in hand we are now confident that the double 
hump is a result of the surface projection of the 
ellipsoidal binary component that fills its Roche 
lobe and orbits its more massive companion on 
a relatively high-inclination orbital plane (to the 
line of sight), coupled with the effect of gravita- 
tional darkening. The difference in minima dips, 
on the other hand, is the consequence of the heat- 
ing of the surface of the Roche lobe-filling star that 
faces a hotter, but much more compact compan- 
ion. This phenomenon is often observed in com- 
pact binaries. It is also referred to as reflection 
(Wilson & Devinney 1971). In what follows, we 
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Fig. 6. — The radial velocity curve of TS 01 folded 
with the orbital period. The points marked in 
magenta are measurements from Gemini spectra. 
The CFHT points are presented in blue. The lat- 
ter were measured anew with the same parameters 
deduced from the line profile fitting by FITSB2. 
The sin fit to the data presented here is practi- 
cally indistinguishable from RV curve calculated 
by Nightfall (see below). 



qualify by cool or optical the Roche lobe-filling 
component, since it is the main contributor of light 
in the optical range, and by hot or X the hotter 
component which irradiates its cooler companion. 
We refrain from the usual wording of primary and 
secondary components in this paper, since, as we 
later discuss, the roles of primary and secondary 
changed during the evolution of this system. 

The optical spectra are too sparse to determine 
the orbital period independently, but they cover 
almost all orbital phases. The orbital period esti- 
mate from the photometry is good enough to cal- 
culate the orbital phases for the spectra. There- 
fore, by measuring the radial velocities of the ab- 
sorption features in each spectrum, we are able to 
construct the radial velocity curve. We used the 
FITSB2 procedure to measure the RV of absorp- 
tion features. FITSB2 performs a simultaneous 
fit of the spectra covering different orbital phases, 
i.e., all available information is combined into the 
parameter determination procedure (Napiwotzki 
et al. 2004). The fit results are the stellar param- 
eters as well as the RVs (Napiwotzki et al. 2004). 
We followed the same procedure as Tovmassian 
et al. (2004) adding the new Gemini spectra. How- 
ever, here, we sought solutions with relatively low 
temperatures, because we now had a better idea of 
the temperature composition of the binary (Tov- 
massian et al. 2007). The best fits were achieved 
with T off = 60 000 ±5 000 K, \ogg = 5. 17 ±0.07 

The radial velocity (RV) curve presented in 
Fig. 6 is fitted with a simple sinusoid. The phase 
zero corresponding to the -/+ crossing of the RV 
curve coincides with the deeper minimum in the 
light curve. It confirms our interpretation that, at 
phase 0.0, the optical component is seen in con- 
junction from behind, turning to the observer its 
smallest projected area and coolest surface tem- 
perature. At phase 0.5, the optical component is 
seen with the same projected area, but presents 
the face with the highest temperature, as result 
of the intense heating from the X component. A 
formal orbital solution leads to semi-amplitude of 
RV K coo i = 216 ± 10 kms -1 , and systematic ve- 
locity of the system 7 — ± 12 kms -1 relative to 
the nebular emission lines. 

The correct determination of phases and RVs 
allowed us to combine all 12 Gemini spectra, elim- 
inating the nebular emission line components, and 
delineating the absorption profiles with increased 
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signal-to-noise. The co-added spectrum is shown 
in Fig. 1. In combination with the UV stellar spec- 
tra from the HST observations, where the stellar 
component is easily separated from the nebular 
one due to the high spatial resolution, the coad- 
ded spectrum allows us to perform a model atmo- 
sphere analysis of the cool component (see Sect. 
3.2. 

3.2. Temperatures and Radii. 

TS 01 is clearly the first planetary nebula 
known to contain a double-degenerate binary (De 
Marco et al. 2008). The first indication that the 
previous interpretations (Tovmassian ct al. 2004; 
Pequignot & Tsamis 2005) of the central star of 
TS 01 may be incorrect came from the shape of its 
light curve (Napiwotzki et al. 2005). Now, armed 
with the X-ray data, we know that even a 130 
kK star can not produce the observed X-ray flux 
and that an additional component is required to 
explain the observed spectral energy distribution 
(SED). 

We analyze the SED, fitting it with a black- 
body as a first approximation. The actual atmo- 
sphere can be significantly different from a black 
body at wavelengths shorter than 900 A, but as 
starting point a black body gives us a good idea 
of what we are dealing with. We will discuss de- 
viations from blackbody later in the paper. The 
data were de-redden according to Schlegel et al. 
(1998) with a canonical ratio of total-to-selective 
absorption, E(B — V) = 0.03 mag and Ry = 3.1. 
We simultaneously fit two black bodies to the ob- 
served spectral energy distribution by introducing 
best guesses for T coo i, T ho t, r coo \/D, and r hot /D, 
where D is the distance to the object, and calcu- 
lating their best fit values. Since there is a gap be- 
tween extreme UV and X-ray wavelength ranges, 
and since the slope of the X-ray data is not strictly 
that of a black body, we obtain three distinct solu- 
tions with similar \ 2 by varying the input param- 
eters. Possible solutions are presented in Table 1. 

The resulting fits are presented in Fig. 7. The 
black body solution for the hot component is 
T hot « 180, 000 ± 25 000 K. Minimum X 2 ~ 2.6 
can be achieved with significantly different tem- 
peratures for the hot component, depending upon 
which part of X-ray data the fit relies on. But 
the hot solutions with T hot > 175 000 K, marked 
in the Fig. 7 as shaded area, do not work, as can 



be seen below, because only certain ratios between 
the hot and cool components fluxes can produce 
the required difference in the depths of minima of 
the light curve. To restrict the range of possible 
solutions, we analyze the form of the light curve, 
together with the radial velocity curve, using the 
binary star modeling code Nightfall 4 . 

Nightfall is based on a physical model that takes 
into account the non-spherical shape of stars in 
close binary systems, as well as the mutual ir- 
radiation of both stars, and a number of addi- 
tional physical effects such as gravitational dark- 
ening and albedo. We fitted simultaneously the 
light curves in three bands and the radial velocity 
curve. The program uses differential magnitudes 
and is tailored to the Johnson photometric system. 
Taking into account that the shape and range of 
amplitudes of the light curve is (a) similar in the 
BUSCA narrow filters and the Johnson V filter 
and (b) docs not show large wavelength depen- 
dence in the optical range, we assigned BUSCA uv 
to the U filter, b to the B filter, and r to the R fil- 
ter. We did not use nir band data, since it was the 
noisiest and would not add anything substantial to 
the analysis. We searched for solutions by setting 
the temperatures to the values estimated from the 
SED. We also fixed the total mass of the system 
close to the Chandrasekhar limit of 1.39 M Q . The 
real M tot might be slightly lower or higher, that 
would not affect this analysis. Leaving the total 
mass parameter free, Nightfall tends to solutions 
involving massive stars, which are excluded. How- 
ever, limiting the total mass still results in a vari- 
ety of remaining parameters that achieve similarly 
good fits. The inclination angle of orbital plane 
to the line of sight, i, must be relatively high to 
produce the observed light variation due to the el- 
lipsoidal form of the cool component, but not too 
high to produce eclipses, which we do not observe, 
neither in the optical, nor in X-rays. Inclination 
angles ranging from the high 40's to the low 70 
degrees are acceptable, and the light curve form 
depends weakly on i within that range. The essen- 
tial parameters for which we seek solutions are the 
masses and radii of the binary components. How- 
ever, the RV and light curves fitting does not pro- 
vide any clue on masses, so additional constrains 



4 http: / / www.hs.uni-hamburg.de/DE /Ins/Pcr/Wichmann / 
Nightfall.html. 
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Fig. 7. — The spectral energy distribution (SED) of TS01. The observed spectra are presented by crosses 
with error bars (SDSS + HST/STIS + FUSE+XMM). The observations are fitted with two blackbodies. 
Three pairs of blackbody solutions are presented according to Table 1. Dashed, dashed-dotted and dotted 
lines represent the cool, the hot components and their sum respectively. Shaded is the area in which the the 
radius of hot component determined from the light curve analysis starts to deviate from the one required for 
SED fitting. In this area the luminosity of hot component does not provide enough irradiation to the cool 
component to produce observed light curve. Thus, only a cool ~ 160 000 K solution is viable. 



9 



Tabic 1: Output parameters of double black-body fits to the SED. 







optical component 






X component 








Solution 


T 


R/D 




T 


R/D 


R a 


log L a 


x 2 




K 


xl0~ 13 


kpc 


K 


xl0~ 13 


R 


L 




Cool 


47 700 


3.80 


21.7 


152 600 


1.20 


0.12 


3.81 


2.54 


Intermediate 


58 900 


3.83 


21.6 


174600 


0.43 


0.04 


3.14 


2.63 


Hot 


60 600 


3.83 


21.5 


205 200 


0.15 


0.014 


2.51 


2.65 



"The values of these columns were calculated assuming R coo i = 0.43R©, and corrected by a factor D ~ \fF = 0.85 because the 
model atmosphere flux F a t m ~ 0.73 X i<BB with T = T e fj at the optical wavelengths, as it follows from the analysis below (see 
3.2.1 and 3.2.2). 



were necessary. 

3.2.1. The cool component. 




5.2 5.1 5 4.9 4.8 4.7 4.6 

log T eff , K 

Fig. 8. — T c ff- log g dependence for post-AGB 
models (Weiss & Ferguson 2009, and priv. 
comm.) Solid line - Z=0.004 for 0.529 & 0.533 M . 
Dashed line - Z=0.0005; 0.539, 0.551 M core- 
spondingly. Tracks for more massive stars spread 
out toward hotter temperatures. The ellipse indi- 
cates range of possible solutions for TS01 stem- 
ming from variety of methods used in this paper. 

The first estimate of the mass of cool com- 
ponent (Tovmassian et al. 2008, 2007) based on 
T c ff and evolutionary tracks for solar composi- 
tion post-AGB stars from (Schoenberner 1983; 
Bloecker 1995) led to M cool w 0.58 M . Recent 
models for different metallicities (Weiss & Fergu- 
son 2009) suggest 0.52 M as the lower limit of the 
mass of a star that heats up to <~ 60 000 K (see 
Fig. 8). For Z=0.001, the mass limit is actually 
0.54 M . Fixing the mass of the cool component 



in Nightfall to that value, we find that regardless 
of other poorly-constrained parameters of the hot 
component, the cool star must have a radius of 
at least 0.42 R to fill its Roche lobe up to 94- 
99% in order to produce the observed light curve. 
Since the cool component is ellipsoidal in shape, 
this radius, as determined by Nightfall, represents 
the mean radius. In fact, the radius depends only 
weakly on the mass adopted for the star, and is 
a stronger function of the binary system's mass 
ratio, which determines the size of Roche lobe. 
For range of mass ratios stemming from the to- 
tal mass of the system in 1.3-1.45 M interval 
the mean radius of the cool component lies within 
0.42-0.45 R range. 

The mass and radius obtained for the cool com- 
ponent leads to a mean log g of 5.03 ±0.03, a value 
deduced by averaging the unevenly distributed 
gravitational acceleration on its surface (Djurase- 
vic 1992). A very similar value of surface gravity 
is obtained by using FITSB2 (§3.1). 

Next, we modeled the stellar atmosphere of the 
cool component using the Tubingen NLTE Model- 
Atmosphere Package (TMAP) 5 (Werner et al. 
2003; Rauch & Deetjen 2003). This code com- 
putes plane-parallel or spherical Non-LTE model 
atmospheres in radiative and hydrostatic equilib- 
rium and considers opacities of all species from 
hydrogen to nickel. The determination of T coo i 
is based on the evaluation of the ionization equi- 
librium through the C in A 1175 A / C iv A 1169 A 
stellar absorption line ratio. We find Tcooi = 
55 ± 5kK. At such a temperature, a surface grav- 



5 http://astro.uni-tuebingen.de/~rauch/TMAP/TMAP.html 
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Fig. 9. — Comparison of a TMAP SED with observation, a. the STIS NUV range. The Hen Fowler series 
(marked at top) are well reproduced; "is" denotes interstellar lines, b. the GEMINI North optical range. 
The H and He lines are marked at top. 



ity of log <7cooi = 4.9 ± 0.5 gives a good agreement 
with the observed spectral line profiles (see Figs. 9 
a & b, comparing a TMAP model with the STIS 
NUV observations and with the gemini spectrum, 
respectively). T coo i and log<? CO oi cannot be better 
constrained given the quality of the data. How- 
ever, the agreement between the values derived 
with TMAP and those derived with FITSB2 and 
Nightfall confirms our correct assessment of the 
basic parameters of the cool component. 

We performed some TMAP test calculations in 
order to derive upper abundance limits for some 
metals. These are summarized in Tab. 2. Note 
that the resonance lines of C IV and N V were not 
used in our abundance determination, since they 
were found to be affected by interstellar line ab- 
sorption. The chemical composition of the cool 
stellar component is close to that of the nebula 
(Stasihska et al. 2009). 

In Fig. 10, the model is compared with observa- 
tions, covering continuously the whole range from 
900 to 10 000 A. Apart from the observations and 
atmospheric model, the nebular emission is shown 
as deduced in Stasinska et al. (2009), and the black 
body curves as implemented in Nightfall. The 
black body corresponding to the cool and hot com- 
ponents are denoted by open stars (red and blue 



Table 2: Photospheric abundances of the cool com- 
ponent of TS01. Ca-Ni are represented by a 
generic model atom. The errors for H, He, and 
C are about 0.3 dex. For N, O, Si, and the iron- 
group elements upper limits are given. 



element mass fraction 



H 




7.471E-01 


He 




2.525E-01 


C 




1.335E-04 


N 


< 


8.306E-05 


O 


< 


7.116E-05 


Si 


< 


6.737E-05 


Ca-Ni 


< 


1.319E-06 
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Fig. 10. — The spectrum of TS 01 from the far UV to the near IR. The observational data after extinc- 
tion correction and nebular emission lines removal are denoted as follows: FUSE - magenta; HST (FUV 
+NUV+UVopt) - blue; SDSS - green. The models for the nebular and stellar emission are shown by grey 
dotted lines. The thick black line is the sum of the model continuum emission from the nebular and stellar 
components. The black-body solutions are presented by red open stars for the cool component and by blue 
stars for the hot component. The thick red line is sum of hot (162 kK) and cool (~ 57 kK) blackbodies plus 
nebular emission. 
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Fig. 11. — From left to right, the Nightfall fits to the light curves in uv, b, and r. In the main panel, the 
red, green, and blue curves represent the solution for X-component with T g of 162 kK, 182 kK, and 202 kK, 
respectively. The deviations with respect to the coolest and hottest solutions are shown in top and bottom 
panels, respectively. The complete set of parameters used in fitting is given in Table 3. The simultaneous fit 
to the RV curve is similar to the sin curve presented in Fig. 6. 



respectively). The sum of the two black bodies 
and the nebular continuum within the observing 
slit as computed in Stasihska et al. (2009) is rep- 
resented with the thick red line in Fig. 10. The fit 
of the models to the observations is excellent from 
the near infrared to 1500 A. 

The curve representing the stellar model plus 
the nebular emission departs slightly from the ob- 
servations at shorter wavelengths. The problem 
with the flux above 1500 A in the NLTE mod- 
els has been noted before (Rauch 2008). More 
importantly, the model is calculated for a spher- 
ically symmetric star with a homogeneous tem- 
perature distribution over its surface. However, 
we know that the cool component is strongly ir- 
radiated and gravitationally distorted, which af- 
fect both its SED and the gravitational acceler- 
ation over the surface of the star. We approx- 
imate the observed, phase-averaged spectrum by 
the non-irradiated atmosphere model even though, 
in some orbital phases, we observe the irradiated 
hemisphere of the cool component. The spectra 
of irradiated atmospheres are flatter in the far UV 
in comparison with the spectra of non-irradiated 
atmospheres with the same T e g (see next subsec- 
tion). Therefore, the average flux in the far UV 
decreases in comparison of the non-irradiated at- 
mosphere spectrum. Overall, we have very good 



agreement between models and observations. 

3.2.2. The hot component. 

The parameters of the hot component are less 
certain. Our knowledge of the hot component is 
based on the binary period, the fact that the cool 
component is partially irradiated to produce the 
observed light curve, and the X-ray flux, which 
cannot originate from the cool component. There 
is an extensive argument in Napiwotzki et al. 
(2005) discussing why the hot component should 
be a compact object with a mass exceeding the 
mass of the cool component. However, neither the 
light curve, nor the X-ray spectra allow us to de- 
termine the temperature or radius of the hot com- 
ponent as well as we did for the cool component. 

We face two problems in the case of TS 01 's hot 
component. First, TS01 is only detected in short, 
soft end of the X-ray range. Second, the calibra- 
tion of data at the extreme soft end of the XMM- 
Newton spectral range is not very reliable. Both 
prevent us from fitting exact atmospheric models 
to it. Hence, like most studies of supersoft X- 
ray sources, we are forced to continue the analy- 
sis using the blackbody that successfully describes 
the spectrum of the hot component in the optical 
and UV range up to 4 x 10 15 Hz (see Fig. 7). We 
also note that this introduces an overestimate of 
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Tabic 3: Parameters of Nightfall fits. 







opt/UV 


component 








X 


component 






Solution 


T 


M 


R 


log g 


inc 


T 


M 


R 


log g 


x 2 




K 


M 


R 




deg 


K 


M Q 


Rq 






Cool 


57100 


0.537 


0.44 


5.02 


52.8 


162 195 


0.853 


0.135 


6.1 


46.27 


Intermediate 


57264 


0.537 


0.43 


5.04 


53.0 


182 195 


0.853 


0.113 


6.3 


46.37 


Hot 


57104 


0.537 


0.43 


5.04 


53.2 


202 000 


0.853 


0.088 


6.5 


46.43 



the luminosity of the X-ray source (Heise et al. 
1994; Swartz et al. 2002). On the other hand the 
temperature can be either overestimated (Swartz 
et al. 2002), or underestimated (Heise et al. 1994; 
Ibragimov et al. 2003). Therefore, when estimat- 
ing the temperature of the hot component in the 
optical/UV range separately from estimates in X- 
rays, we should not worry too much if discrepan- 
cies arise. It is within the optical/UV range that 
the irradiation of the cool component matters, so 
we may fix the parameters of the cool component 
in Nightfall and seek solutions for the hot com- 
ponent. Even so, varying freely both the radius 
and temperature of hot component, we still do not 
reach unambiguous solutions. 

In Fig. 11, we present the fits produced by 
Nightfall to light curves in three different filters. 
The parameters for the fits for the different tem- 
peratures of the hot component are presented in 
Table 3. The temperature and the mass of the cool 
component were kept fixed (the difference in mean 
temperatures of the cool component in the Tabic 
reflects the different degree of irradiation). The 
parameters that were fitted are the orbital inclina- 
tion angle, the fraction to which both components 
fill their Roche lobes, and the temperature of the 
hot component. Three models with different tem- 
peratures for the hot component are displayed in 
Fig. 11. The fits shown red, green, and blue corre- 
spond to Thot= 162, 182, and 202 kK, respectively. 
The lower and upper panels show the deviations of 
the fit from the observations for the extreme cold 
and hot solutions. As can be seen from Fig. 11 
and Table 3 the differences between three models 
ranging from 160 to 200 kK are not important in 
the optical domain. 

Note, that the range of temperatures for hot 
components obtained from blackbody fitting to 
SED and from the fitting of light curves by Night- 
fall are similar. But comparison of Tables 1 and 
3 shows that only for low temperature solution 



the radii deduced by both methods are compat- 
ible. Introducing the radius of the cool compo- 
nent deduced from Nightfall into the R coo i/D pa- 
rameter used in the fit of two black bodies to 
the SED results in a distance of 25 kpc for the 
minimum radius of R coo ; = 0.42 R Q and leads 
to i?hot ~ O.13i?0 for a blackbody temperature 
150 kK. The color temperature of the hot compo- 
nent is probably slightly higher. This is caused 
by the divergence of the real atmosphere from the 
black-body at high energies and also by the use 
of a Roche lobe-shaped cool component in Night- 
fall, instead of the spherical shape in all other cal- 
culations. Similarly, using stellar atmospheres in- 
stead of black-bodies reduces the distance to about 
21 kpc, as reflected in Table 1. Increasing tem- 
perature is compensated naturally by a smaller 
radius. The solutions with temperatures above 
175 kK come up with parameter Rhot/D too small 
to be compatible with results of light curve fitting. 
In Fig. 7 the blackbody solutions of hot component 
that fall into shaded area are not luminous enough 
to provide necessary irradiation and produce the 
observed light curves. In the meantime, tempera- 
tures above 185 kK are not tolerated by ionization 
modeling of the nebula (Stasihska et al. 2009). 

Therefore, we consider 160-175 kK temperature 
range and Rhot ~ O.IR© to be the closest to the 
real properties of the hot component. Note that 
the hot component has a radius Rhot > O.O4R 
at least. This is much larger than an ordinary 
white dwarf. As such TSOl's hot component is 
very similar to the supersoft X-ray source (SSS) 
Lin 358, one of the two SMC symbiotic stars stud- 
ied by Orio et al. (2007). Majority of estimates of 
temperatures and luminosities of super-soft X-ray 
sources are made using black bodies and values 
obtained here are useful when comparing to other 
similar objects. 

Fitting an exact atmospheric model to the X- 
ray data for the hot component does not make 
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Fig. 12. — Comparison of the observed light curves 
(filled symbols) at A1350A (top panel), A1650A 
(middle panel), and A1900A (bottom panel) with 
the computed light curves in the same bands (solid 
and dashed lines) . The horizontal error bars mark 
exposure length for each point. The light curves at 
A1650A are obtained separately from HST FUV 
(circles) and NUV (triangles) detectors. Their 
disparity indicates the difference of light curve 
smoothing due to the exposure length. The model 
fluxes, integrated in the same phase ranges as 
observed ones are shown by open symbols. All 
models are computed for the cool solution from 
Table 2. To account for differences between the 
black body and real atmosphere we introduced 
the parameters a & b (see text), obtained empiri- 
cally from comparison of these two observed light 
curves. The values of parameters for flux correction 
in corresponding bands are marked in the plot. 

much sense, because the observed energy range is 
too small and the quality of the data is too poor. 6 



For the analysis of the ionization of the nebula, however, 
using a black body would cause serious problems at high 
energies. This is why, for their photoionization modelling, 
Stasinska et al. (2009) selected a suitable spectrum from 



An additional test for checking the estimated 
temperature of the hot component comes from 
modeling of the light curve in the UV. However, 
Nightfall cannot calculate model fluxes in the UV, 
so we implemented our own code (Shimansky et al. 
2002), which calculates the irradiation in a binary 
according to the prescription in Howarth & Wil- 
son (1983). Nightfall is based on the same algo- 
rithm, so we naturally obtained exactly the same 
light curves for optical bands. We compute the 
UV light curves for the cool soloution from Table 
3, and present it, together with the observed one, 
in Fig. 12. In computing the UV light curves we 
take into account the difference between a black 
body and a stellar atmosphere introducing the ra- 
tio f x = F A (BB)/F A (BSA) of black body and a 
stellar atmosphere fluxes at a given wavelength 7 . 

3.3. The total mass. 

The mass of the hot component cannot be de- 
termined directly from the observed data. How- 
ever, simultaneously fitting the RV curve and light 
curves indicates a tendency towards improvement 
as the total mass approaches the Chandrasekhar 
limit. In Fig. 13, we plot the calculated \ 2 as a 
function of the total mass of the system, by fixing 
in Nightfall the temperatures of the components 
to our best estimates, i.e. 57 and 162 kK and 
adopting a mass of 0.54M Q for the cool compo- 
nent. The lower limit for total mass is of 1.25 M Q , 
corresponding to a hot component with a mass of 
~ 0.7 M . This mass corresponds to the lower 
limit for a white dwarf to sustain steady nuclear 



a grid of models with halo composition, which reasonably 
well describes the observed fluxes both in opt/UV and X- 
ray range and provides a more realistic picture. 
7 In the optical part of the spectrum / opt ~ 0.73 for both 
the cool and hot components; near the Lyman edge the 
SED of the cool component is different from a blackbody: 
/l350 ~ 1-18 and /ibso ~ 1-02; the values of fx for the hot 
component remain close to 0.73. Moreover, the values of 
fx for the cooJ component must depend on the irradiation 
flux. If the irradiation flux increases, fx decreases, because 
the spectra of the irradiated stellar atmosphere are getting 
closer to a black body spectrum. The best approximation of 
the observed light curves by the model light curves was ob- 
tained for a simple linear dependence fx = a — b- (Fi rr /Fo), 
where fj rr is the irradiation flux at a given point of the cool 
component surface, and Fo is the flux from the cool com- 
ponent at the same point. A change of the continuum slope 
in the UV band at the various orbital phases, observed by 
HST, confirms this picture. Parameters a and b deduced 
for each A are presented in Fig. 12. 
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Fig. 13. — The line with dots depicts the depen- 
dence of the x 2 value on the total mass of the 
system. Fixed parameters in the calculations were 
the temperatures of the components (55 & 160 kK 
respectively) and M coo / (0.54 M©). The free pa- 
rameters were the inclination of the binary orbit i, 
which changed similar to the x 2 m a range of val- 
ues shown on the upper axes; and the Roche lobe 
filling factors (//) of both components shown as 
shaded areas in a range of values marked on the 
right side axes. The dark shaded area corresponds 
to Roche lobe filling factor of the hot component 
and the light grey to that of the cool component. 



ting the light and RV curves at the same time, but 
of course a lower than Chandrasekhar limit mass 
is not excluded. Better measurements of radial 
velocities are required for a more reliable determi- 
nation of stellar masses in TS01. 
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Fig. 14. — Position of the hot component in the 
temperature-luminosity diagram (filled box) re- 
spective to the tracks of hot accreting white dwarfs 
in the steady-burning approximation (Iben 1982). 
The bold solid curve shows the high-temperature 
boundary of the stable-burning strip. The num- 
bers by the tracks are corresponding white dwarf 
masses in M©. 



burning on its surface. Increasing the total mass 
of the binary from that minimum value produces a 
decrease in the x 2 value, with a small plateau of \ 2 
at about 1.4M©. The % 2 keeps falling as the total 
mass increases. This is due to the fact that ris- 
ing mass of the hot component shrinks the Roche 
lobe of the cool star. At around M to tai = 1-47 the 
cool component fills its corresponding Roche lobe 
to 99.9%, which helps a better fitting of the light 
curves. Improvement of x 2 from there on is con- 
ditioned by the rapid increase of size (Roche lobe 
filling factor = //) of the hot component, which 
might be unrealistic. The mass could in princi- 
ple be constrained by the fit to the RV curve, but 
unfortunately the RV data is too poor in quan- 
tity and quality to have strong influence on x 2 - 
We consider the flattening of the x 2 — Mtotai curve 
around M to tai = 1-40 as an indication of the best 
solution, where a balance is achieved between fit- 



The mass estimate of the hot component can 
be checked from its position in the H-R diagram 
(Suleimanov & Ibragimov 2003). Figure 14 shows 
position of the hot component of TS01 in the log 
T ff — logL diagram. It suggests that the ob- 
ject has a mass of about 0.78 ± 0.02 M , if it 
has not cooled much since active intense accre- 
tion ceased. The lower luminosity boundary of 
the rectangle corresponds to the upper luminos- 
ity limit obtained from the SED, while the upper 
luminosity boundary corresponds to the luminos- 
ity, obtained from the optical light curves model- 
ing. Given all uncertainties, the agreement with 
mass estimate obtained above seems to be rather 
good. But we note cautionary that the tracks were 
calculated assuming solar chemical composition of 
accreted matter, while real chemical composition 
of the hot component is not known and that the 
structure of its envelope, which defines the rate of 
cooling, may be different from the one of a steady- 
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burning object. 

For the readers convenience we summarize pa- 
rameters of the binary estimated by variety of 
methods in Table 4. 

4. Evolutionary considerations. 
4.1. Formation of TS 01 

Below, we present examples of evolutionary sce- 
narios that can result in the formation of a sys- 
tem reasonably similar to TS01. These scenarios 
envision the formation of a massive white dwarf 
that cools for a long time but is reheated by com- 
pression due to accretion and nuclear burning of 
material captured from the stellar wind of com- 
panion. The core of the latter is currently the less 
massive component of the binary nucleus of TS 01. 
In the analysis of the origin of TS01, one has to 
take into account its location in the Galactic halo, 
where star formation ceased about 10 Gyr ago (e. 
g. Marin- Franch et al. 2009). This sets an upper 
limit of < 1 M for the mass of the progenitor of 
the cool component. 

The presence of an old white dwarf and a 
nascent white dwarf in a system with an orbital 
period of 3.92 hr only implies that previous evolu- 
tion of the system involved common cnvelope(s). 
First, let us assume that both components of the 
nucleus of TS01 had AGB precursors. This sug- 
gests the following scenario 8 . The initial masses 
of the components are significantly different. The 
initial system is so wide that the primary evolves 
to the AGB unaffected by the presence of the com- 
panion. In the FGB 9 and AGB the system might 
manifest itself symbiotic system with AGB 
and main-sequence components (Kenyon & Web- 
bink 1984). Owing to the wind mass loss from 
the system, the separation of the components in- 
creases. Both the radius of the primary and the 
radius of its Roche lobe increase, with the radius of 
the primary growing faster until the primary over- 
flows its Roche lobe (RLOF) close to the tip of the 
AGB. Both because the primary at this time has 



Below, we use the terms "primary" and "secondary" for 
components that were, respectively, more and less massive 
at ZAMS. As we will show further, "primary" is precursor 
of the hot component, while "secondary" is precursor of the 
cool component. 

The first giant branch, in which the helium nucleus is 
formed but not burning yet. 



a deep convective envelope and the mass ratio of 
the components is high, the dynamical mass loss 
is unavoidable (Hjellming & Webbink 1987) and 
the shedding of the envelope results in formation 
of a common envelope (CE) and a reduction of 
the separation of the components due to angular 
momentum loss in CE. What remains of the pri- 
mary after this episode is the more massive com- 
ponent of the core of TS01. The system remains 
wide enough so that other component may evolve 
to become a giant star too and experience RLOF 
close to the tip of the AGB, forming the current 
cool component. The matter ejected during the 
second CE episode is now observed as a planetary 
nebula. 

We now present numerical estimates that ar- 
gue in favor of the feasibility of a scenario such as 
that just described. In our evolutionary simula- 
tions, we use the "rapid evolutionary code" SSE 
(Hurley ct al. 2000) based on the analytical fits 
to detailed grids of full stellar models. We use 
the SSE code because detailed evolutionary tracks 
for low-metallicity stars with M < 1 M have yet 
to be computed. Comparison with data for more 
massive stars (e. g. Weiss & Ferguson 2009) 
shows that the initial-final mass relations used by 
us agree with the results of sophisticated, full evo- 
lutionary computations to within ~ 10 per cent 
and, hence, qualitatively, the resulting scenario 
must be robust. We note also, that the results 
of evolutionary calculations depend heavily on the 
opacities used for the models. 

Stasihska ct al. (2009) estimate that the metal- 
licity of TS01 ranges from 1/12 to 1/30 of the 
solar value (taken as Z =O.O14; Lodders et al. 
2009). For our calculations we accepted Z=0.001 
as a proxy to the metallicity of TS01, since our 
goal is to demonstrate the possibility of forming 
a system similar to TS 01, rather than attempt to 
reproduce precise values for parameters that are 
still quite uncertain. 

The second CE episode, which produced the 
current cool component, followed RLOF by its 
precursor close to the tip of AGB. Using SSE, we 
find that, for Z=0.001, stars with a ZAMS mass 
exceeding 0.89 M evolve to the tip of AGB in less 
than 10 Gyr. At the tip of the AGB, a star with 
M Z ams = 0.89 M has a mass of 0.60 M and 
a core mass of 0.54 M , which is coincidentally 
similar to the estimated mass of the cool compo- 
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Table 4: Summary of parameters of TS01. 

T C ooi M coo i R coo i logg coo i T hot M hot R hot used data 

Method xl0e3K M R xl0e3K M R 

FITSB2 601T5 5.17 ±0.07 GEMINI, CFHT 

SED 2BB fit 55 ± 7 0.43* 152-175 0.04-0.12 SDSS, HST, FUSE, XMM 

Evolution 0.53 ± 0.01 Weiss k Ferguson (2009) 

NIGHTFALL 57 * 0.54 * 0.43 ± 0.02 5.03 ±0.02 160-200 0.71-0.93 0.1-0.13 Calar-Alto, GEMINI, CFHT 

TMAP 55 ±5 4.9 ±0.5 GEMINI, HST 



* Fixed parameters in corresponding fits 



nent in TS01. Motivated by Fig. 13, we adopt a 
total system mass of 1.39 M . Given a mass of 
0.54 M for the the cool component, the mass of 
hot component is then s=y 0.85 M Q , corresponding 
to an initial mass of 2.5 M . 

At the tip of AGB, the precursor of the cool 
component had a radius R ~ 140 R Q . Using the 
formula from Eggleton (1983) for the dimension- 
less radius of the Roche lobe r L , we estimate that 
the pre-CE separation of components was about 
400 R . The variation of the separation of com- 
ponents in CEs may be described by the formula 
suggested by Webbink (1984): 

a f M r L 2 M 2 - M r n ~ 



Mc 
M 2 



1 + 



2.L 



M 1 



(1) 



where M 2 and M c are initial and final masses 
of mass- losing component (the donor), Mi is the 
mass of companion, a ce \ is the product of effi- 
ciency of common envelope expulsion a ce and the 
structural parameter A which characterizes bind- 
ing energy of the donor envelope, r 2 .L is the frac- 
tional Roche lobe radius of the donor. The reduc- 
tion of the separation from an. ~ 400 R to the 
current a/ ~ 1.3 R© is possible if a ce X w 0.0015, 
i.e., is extremely low. 

In the stage preceding the common envelope, 
the system contained an AGB star and a massive 
(hot) companion accreting from the wind. In this 
stage, the system could be identified with a sym- 
biotic star (Tutukov & Yungelson 1976; Kenyon & 
Webbink 1984; Yungelson ct al. 1995; Lu et al. 
2006). Accretion reheated white dwarf and re- 
sulted initially in unstable and later in stable hy- 
drogen burning at the surface of white dwarf. En- 
ergy release by nuclear burning also contributed 
to the heating of white dwarf. 

It is plausible that currently hot component 
still burns remainders of hydrogen accreted in this 



stage. During the symbiotic stage, the precursor 
of the cool component lost about 0.29 M Q via a 
wind. Accretion from the wind in symbiotic sys- 
tems is inefficient (~ 10% de Val-Borro ct al. 2009) 
and we may safely assume that all mass lost by 
the donor was lost from the system taking away 
specific angular momentum of the donor ( "Jeans 
mode of mass ejection") and that the mass of the 
hot component did not change. Jeans mode of 
mass ejection has an invariant a x (Mi + M 2 ) and, 
hence, the separation of the components in the be- 
ginning of the symbiotic stage was 320 R . This 
separation, 320 R , is also the separation of com- 
ponents after the first CE stage, which aborted 
the ascend of AGB by the initially more massive 
component close to the tip of the AGB. Before the 
first CE stage, the mass of the star decreased via 
wind mass-loss from 2.5 M to 1.29 M©. Its ra- 
dius at the tip of the AGB was 495 R . Like for 
the second RLOF episode, from the condition of 
RLOF we may estimate that the separation of the 
components at the beginning of the first RLOF 
was 1200 R©. The reduction of the separation in 
the first CE phase from 1200 R to 320 R implies 
a ce X ~ 1.5. The first CE episode could also have 
been preceded by a symbiotic stage. We again as- 
sume that all mass lost by the donor was lost from 
the system via the Jeans mode of mass ejection. 
We then estimate the initial separation of com- 
ponents as close to 770 R©. We neglect the wind 
mass loss during the first red giant stage, which is 
only several 0.01 M© for M ~ 2.5 M . The nu- 
merical data is summarized in the upper part of 
Table 5 (scenario I) and is presented in a form of 
a cartoon in Fig. 15. 

An apparent problem with the suggested sce- 
nario is the large difference in a ce X for the com- 
mon envelope stages. While a ce \ <~ 1 is typical for 
WD+MS stars with periods below about 10 days, 
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Table 5: Numerical evolutionary scenarios for TS01. The upper part of the Table shows scenario (I) based 
on energy balance formalism of Webbink (1984), while lower part of the Table presents scenario (II) based 
on angular momentum balance formalism of Nclcmans ct al. (2000) 



Mi 


M 2 


a 


Comment 


M 


M 


R© 




2.5 


0.89 


770 


ZAMS 


1.29 


0.89 


1200 


The end of AGB ascend by the primary, beginning of the first RLOF (CE) 


0.86 


0.89 


320 


The end of the first CE, a ce w 1.5, formation of the first WD, 








beginning of the symbiotic stage 


0.86 


0.60 


400 


The end of the symbiotic stage; RLOF (CE); ejection of PN, a ce X ~ 0.001 


0.86 


0.54 


1.3 


Present state 



5.0 


0.89 


150 


ZAMS 


5.0 


0.89 


150 


The end of FGB ascend by the primary, beginning of the first RLOF (CE) 


0.87 


0.89 


240 


The end of the first CE, 7 w 1.2, 








formation of He-star evolving into first WD, beginning of the symbiotic stage 


0.87 


0.73 


260 


The end of secondary evolution in E-AGB; RLOF (CE); ejection of PN, a ce \ ~ 


0.87 


0.53 


1.3 


Present state 



which are supposed to form via one common enve- 
lope stage and implies that the energy spent on the 
expulsion of the common envelope is comparable 
to the orbital energy of the initial binary (Nclc- 
mans & Tout 2005), a ce \ ~ 0.001 during the sec- 
ond CE episode appears atypically low. However, 
common envelopes remain virtually terra incognita 
in stellar evolution and we cannot exclude a signif- 
icant difference in the interaction of the AGB star 
envelope with a MS or a WD companion, which 
differ in structure and, most importantly, by two 
orders of magnitude in radius (whereas the drag 
force is <x R 2 ). 

An alternative scenario for TS 01 assumes that 
present cool component had a precursor with 
ZAMS mass of 0.89 M , while the hot compo- 
nent descended from a helium star which was 
formed by RLOF close to the tip of FGB. For 
instance, a 5 M star has a maximum He-core 
mass of 0.87 M which, presumably, evolves into 
a CO WD of the same mass 10 . Thus, the ini- 



'We set the mass of WD equal to the mass of its He-star 
precursor, thus implicitly neglecting the possibility of re- 
expansion of the He-star after exhaustion of helium in its 
core. Such an expansion with formation of a shallow CE 



tial system could contain a 5M component and 
a 0.89 M component and after the 1st CE to be- 
come a (0.87+0.89) M system. If 5 M star filled 
Roche lobe close to the tip of FGB when it radius 
was close to 80 R , prior to RLOF separation of 
components have had to be close to 150 R . 

If the precursor of the cool component was an 
AGB star, the smallest radius with which it could 
overfill its Roche lobe at E-AGB was « 100 R . 
At this moment the total mass of the star was 

0. 73 M©, the mass of the core - 0.53 M©, and from 
the RLOF condition we obtain that the separa- 
tion of the stars was w 260 R . In the second CE 
the separation decreased from 260 R to 1.5 R by 
ejection of 0.2 M . This is possible if a ce X w 0.01, 

1. e. an order of magnitude larger than in the first 
scenario. 

If we account for Jeans-mode mass loss by the 
precursor of the cool component we obtain that, 
after the first CE, the separation of components 
was close to 240 R . Thus we arrive to an appar- 
ent controversy: in the suggested scenario, in the 



and almost negligible mass loss was discovered by Iben & 
Tutukov (1985) for solar metallicity stars, but its possibility 
was never explored for non-solar metallicitics. 
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Fig. 15. — Cartoon depicting the evolution of 
TS01. 



first CE episode, the separation of the components 
had to increase from 150 to 240 R ! 

However, it was noticed by Nelemans et al. 
(2000) and later confirmed by Nelemans & Tout 
(2005) that using Eq. (1) for the description of 
the outcome of unstable mass exchange between a 
giant and a MS-star often does not allow to repro- 
duce well measured parameters of many post-CE 
binaries. As an alternative, Nelemans et al. (2000) 
suggested to estimate the post-CE separations of 
components using an equation for angular momen- 
tum balance: 



Ji - Jf = jJi 



AM 



(2) 



Here J is the orbital angular momentum, sub- 
scripts i and f denote the initial and final values of 
the momentum, AM is the mass lost from the sys- 
tem (the envelope of the donor) , and M tot is total 
initial mass of the system. Thus, a single parame- 
ter 7 describes the fraction of initial specific orbital 
angular momentum of the binary taken away by 
outflowing matter. This "7-formalism" leads to: 



Of 
0(1 



Mi 

M r 



M c + M 2 
Aft,* 



1-7 



AM V 
M^J 



(3) 



Here M c is the mass of the core of the mass-losing 
component. An increase of the separation during 
the CE from 150 R Q to 240 R Q is possible if 7 rj 
1.2. Rather similar combinations of 7 for the first 
CE and a re X for the second one were found for 



some systems studied by Nelemans & Tout (2005, 
see their Figs. 1 and 5). 

4.2. Common Envelope Remnant vs. Sin- 
gle Star Evolution through post- AGB 
phase. 

Ejection of a common envelope definitely dif- 
fers from formation of a planetary nebula by the 
usually assumed superwind mechanism. In that 
context, it is interesting to compare parameters 
of the cool star deduced here with the evolution- 
ary models for post-AGB stars. Given that the 
cool component of TS01 nearly fills its Roche lobe 
and that it is currently contracting, it has only 
recently terminated the phase of common enve- 
lope evolution. The structure and mass of its 
envelope might be very different from that of a 
single star passing through the early epochs of 
planetary nebula nucleus stage. Here we compare 
the derived parameters of the cool component of 
TS 01 with two sets of models of remnants of sin- 
gle stars with initial masses of 1.0 M (lower pro- 
genitor masses are not available in the literature) . 
Based on estimated abundances, we selected the 
models M = 0.623 M , Z=0.001 Vassiliadis & 
Wood (1994) and M = 0.547 M , Z=0.0005 Weiss 
& Ferguson (2009). These models agree well re- 
garding the time-dependence of heating of the core 
of a post-AGB star (lower right panel of Fig. 16). 
TS 01 has T ff similar to them at the age of w 6000 
yr. This age estimate is in a good agreement with 
the age deduced from the expansion velocity and 
distance to the planetary nebula (Stasihska ct al. 
2009). 

Compared to the most modern and the clos- 
est in mass model of Weiss & Ferguson (2009), 
the nucleus of TS 01 is slightly more compact (by 
0.05 dex) and significantly (by more than 0.3 dex) 
less luminous. Since the main source of luminos- 
ity of post-AGB stars is hydrogen burning, this 
may mean that common envelope remnants may 
have less massive H/He envelopes around degen- 
erate cores than their post-AGB counterparts. 

This comparison clearly indicates that evolu- 
tion in common envelopes might alter evolution 
of stars in close binary systems compared to sin- 
gle stars and a more complete analysis of TS01 
is warranted than can be made with single star 
models. 
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Fig. 16. — Various parameters of the models of 
post-AGB stars vs. time. Solid line is for the 
M = 0.547 M Q , Z=0.0005 model Weiss & Fergu- 
son (2009), dotted line is for the M = 0.623 M Q , 
Z=0.001model Vassiliadis & Wood (1994). TS01 
is marked by squares with error bars at age of 
(~ 6000 yr), appropriate to its effective tempera- 
ture. 



4.3. TS01 and SNe la 

The evolutionary path suggested for TS01 in- 
cludes a stage of a symbiotic star which is consid- 
ered as one of the routes to SN la (e. g. Tutukov & 
Yungelson 1976; Iben & Tutukov 1984; Munari & 
Renzini 1992). However, conditions in symbiotic 
systems are not favorable for an efficient accumu- 
lation of matter by the white dwarf components. 
Accretion from the wind typically allows only sev- 
eral per cent of the mass lost by the donor to be ac- 
creted. In the numerical scenario above, the max- 
imum mass-loss rate by the progenitor of the cool 
component estimated by means of SSE is close to 
2 x 10~ 7 M Q yr" 1 . This means that for about 
10 Gyr the white-dwarf (hot) component stays in 
the regime of unstable thermonuclear burning (of 
Novae eruptions) (Nomoto 1982) and instead of 
accumulating mass it may erode. Conditions for 
accretion "improve" if the accretor is located in 
the zone of acceleration of the stellar wind, which 



requires the proximity of the donor surface to the 
Roche lobe (Yungelson et al. 1995), or if the stellar 
wind is pumped close to the Roche lobe by pulsa- 
tions and still remains slow (Podsiadlowski & Mo- 
hamed 2007). Then accretion efficiency may be- 
come close to 100%. Using Bondi & Hoyle (1944) 
formalism for wind accretion, accounting for pos- 
sible location of accretor in the wind acceleration 
zone and taking accretion rate limits for stable 
hydrogen burning after Nomoto (1982), we esti- 
mate that the system could accrete steady for the 
last several 100 000 yr prior to CE and accumulate 
only several 0.01 M Q . In a more general context, 
the circumstances listed above prevent symbiotic 
stars from being efficient progenitors of SN la and 
the estimated rate of occurrence of SN la in these 
systems is only ~ 10 -6 yr" 1 on a Galactic scale (e. 
g. Yungelson 2005). 

For a fraction of time between Novae eruptions 
and in steady-burning regime the system can man- 
ifest itself as a supersoft X-ray source (e. g van 
den Heuvel et al. 1992; Truran & Glasner 1995; 
Yungelson et al. 1996). 

Note that, during the stage of accretion onto 
the current hot component, the matter could in- 
flow onto the equatorial regions of the dwarf while 
it outflows from the polar regions. The "bars" seen 
in TSOl's nebular shell (Stasihska et al. 2009) 
may be the remnants of jets that once existed in 
the system. 

Figure 17 shows the positions of double- 
degenerate systems with known parameters in the 
M tot — P or b plane. As well, positions of sev- 
eral sdB stars with white dwarf companions arc 
shown. The latter systems will turn into double- 
degenerates after completion of helium burning in 
sdB stars. Thus, its short orbital period of 3.92 
hr and its total mass close to the Chandrasekhar 
mass makes TS01 very promising candidate pro- 
genitor for a SNIa in the double-degenerate sce- 
nario for these events (Tutukov & Yungelson 1981; 
Iben & Tutukov 1984; Webbink 1984). For in- 
stance, the merger of components will occur in 
«660 Myr in the first evolutionary scenario sug- 
gested above and in ssl.2 Gyr in the second 
scenario. The only "competitor" to TS01 is 
sdB+WD system KPD 1930+2752 with orbital 
period 2.28 hr, M sdB = 0.45 - 0.52 M , M tot = 
1.36 - 1.48 M Q (Geier et al. 2007). In the lat- 
ter system, subdwarf star will turn into a WD in 
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w(220 - 140) Myr, see Yungelson (2008) for esti- 
mates of lifetime of sdB stars. It will take two WD 
several tens of Myr more to merge. Favourable 
conditions for central carbon ignition may come 
to fruition just in systems with low mass ratios 
of components (Yoon et al. 2007), like TS01 and 
KPD 1930+2752. 
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Fig. 17. — Positions of the known de- 
tached double-degenerate systems (DD) and sub- 
dwarf+white dwarf systems which will turn into 
double degenerates, dl - DD with two visible 
spectra, si - DD with one visible spectrum. Up- 
dated figure from Geier et al. (2007). The position 
of TS 01 is marked by an arrow, indicating the ten- 
dency of better fit to the observed data with in- 
crease of total mass of the binary. The minimum 
value M to t > 1.35M Q corresponds to the leveling 
of x 2 slope in Fig. 13. We consider the values at 
plateau in the \ 2 ~ -^tot dependance as a range of 
best solutions. 



5. Conclusions. 

After a decade of intense study, we have 
achieved a good understanding of an object whose 
discovery spectrum was misidentified and incom- 
prehensible in 1997. Since then, the object has 
been observed at practically all wavelengths with 



the help of the most advanced instruments. This 
paper accompanies Stasihska et al. (2009), which 
focuses upon the chemical composition and ioniza- 
tion state of TSOl's nebular shell. Here, we focus 
on the nature of the close binary nucleus of the 
PN. 

TS01 is one of the shortest period sys- 
tems among the double-degenerate or pre-double- 
degenerate systems, with an orbital period of 3.924 
hours This fact would not have caused confusion 
if the older of the components were significantly 
cooler than the core of the star that most re- 
cently ejected its envelope to form the current PN. 
However, observations and analysis clearly demon- 
strate that TSOl's nucleus is comprised of two 
compact stars, both extremely hot and thus, both 
being sources of ionization for the nebula. This 
unusual phenomenon created confusion and mis- 
interpretation of the object in the past. Neverthe- 
less, the correct understanding of the ionization 
source does not change the essence of those pre- 
vious interpretations. TS 01 remains a PN with 
a record low oxygen abundance (Stasihska et al. 
2009). 

According to our scenario, TS01 evolved 
through two common envelope episodes. In the 
current stage we are observing the remainders of 
the second common envelope as a PN. The core of 
the envelope-shedding post-AGB star is in the pro- 
cess of contraction and heating up. At the present 
time, it nearly fills its Roche lobe and has an el- 
lipsoidal shape. Before the last CE episode, the 
more massive component, which became a white 
dwarf earlier, underwent a period during which it 
accreted mass at a high rate and burned hydrogen 
steady. Since then it stays close to the temper- 
atures range typical for supersoft X-ray sources. 
Its properties make TS 01 one of the softest X-ray 
sources ever, similar to Lin 358 (Orio et al. 2007). 

The parameters of the binary system were de- 
duced using a wealth of information and via three 
independent routes. Although, each of these 
methods requires its own assumptions and each 
alone produces ambiguous results, in combination, 
they converge to values with unusual precision. 
Using the spectral energy distribution, from the 
far infrared to X-rays, the light and radial velocity 
curves, and by fitting atmospheric models to the 
stellar absorption features of the cool component, 
we find that the cool component has a mass of 
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O.54±O.2M , an average T eff of 58 000±3 000 K, a 
mean radius of 0.43±0. 3 R , and logy = 5.0±0.3. 
The cool component nearly fills its Roche lobe. 
The temperature and gravity over the surface of 
the cool component are not homogeneous. 

The chemical composition of the cool compo- 
nent from atmosphere model fitting was deter- 
mined as: 12 + log He/H = 10.95 and 12 + log 
C/H = 7.20, with an uncertainty of about 0.3 dex, 
and upper limits 12 + log N/H < 6.92 and 12 + 
log O/H < 6.80. Overall, the agreement with the 
abundances found in the nebula by Stasihska et al. 
(2009) is very good, except for the carbon abun- 
dance, which is found to be higher in the nebula 
for a reason yet not understood. 

The parameters for the hot component are less 
certain. It is fairly clear that the spectral energy 
distributions of real stars at such high tempera- 
tures depart from that of a black body. The range 
of temperatures that we determined for the hot 
component spans 160-200 kK. It seems that the 
real object acts like a 180-200 kK blackbody in 
the X-ray range but appears as a 160 kK black- 
body in the UV/optical range. Uncertainty in its 
temperature leads to uncertainty in its size, but it 
is obvious from our calculations that the hot com- 
ponent is larger than normal for a white dwarf, 
R-hot > 0.1 R , and is probably bloated as a result 
of intense accretion in the recent past. However, 
we have indirect information on the hot compo- 
nent through photoionization modeling by repro- 
ducing the intensities of the lines emitted by the 
nebula (Stasihska et al. 2009). We estimate the 
distance to the object as <~ 21 kpc, and our most 
reasonable luminosity estimate for the X-ray com- 
ponent is <~ 1O 4 L , appropriate for a supersoft 
X-ray source. 

The total mass of the binary is very close to 
Chandrasckhar limit. This makes TS 01 one of the 
best of the known candidates for the progenitor of 
a type la supernova. 
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